CWRU-P18-02 
astro-ph/02 12083 



Gravitational Leakage into Extra Dimensions: 
Probing Dark Energy Using Local Gravity 

Arthur LueB and Glenn Starkmani 

Department of Physics 
Case Western Reserve University 
Cleveland, OH 44106-7079 



Abstract 



The braneworld model of Dvali-Gabadadze-Porrati (DGP) is a theory where 
gravity is modified at large distances by the arrested leakage of gravitons off 
our four-dimensional universe. Cosmology in this model has been shown to 
support both "conventional" and exotic explanations of the dark energy re- 
sponsible for today's cosmic acceleration. We present new results for the gravi- 
tational field of a clustered matter source on the background of an accelerating 
universe in DGP braneworld gravity, and articulate how these results differ 
from those of general relativity. In particular, we show that orbits nearby a 
mass source suffer a universal anomalous precession as large as ±5 /ias/year, 
dependent only on the graviton's effective linewidth and the global geometry 
of the full, five-dimensional universe. Thus, this theory off'ers a local gravity 
correction sensitive to factors that dictate cosmological history. 
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I. INTRODUCTION 



The gravity theory of Dvah-Gabadadze-Porrati (DGP) is a braneworld theory with a 
metastable four-dimensional graviton The graviton is pinned to a four- dimensional 
braneworld by intrinsic curvature terms induced by quantum matter fluctuations; but as 
it propagates over large distances, the graviton eventually evaporates off the brane into an 
infinite volume, five- dimensional Minkowski bulk. There exists a single free parameter in 
DGP braneworld gravity, the crossover scale, tq. This scale dictates that distance below 
which gravity is controlled by brane effects, but larger than which gravity assumes a five- 
dimensional behavior. As a result, the DGP braneworld theory is a model in a class of 
theories in which gravity deviates from conventional Einstein gravity not at short distances 
(as in more familiar braneworld theories), but rather at long distances. Such a model has 
both intriguing phenomenological as well as cosmological consequences @J|,|,^||,||,|lO 



A braneworld model of the sort where gravity is modified at extremely large scales is mo- 
tivated by the desire to ascertain how our understanding of cosmology may be refined by 
the presence of extra dimensions. Indeed, there exist novel cosmologies in this theory that 
provide an alternative explanation of the cosmic dark energy [Q. 

Recovery of Einstein gravity at short distance scales in DGP [jTT[, especially for static 



sources |ll2| , [l3| , [I^ , p!5| , is a subtle effect. Even though gravity is four-dimensional at distances 



less than tq, it is not always Einstein gravity. For a point source whose Schwarzschild radius 
is Tg, general relativity is only recovered for distances shorter than 

1/3 



{rlr.y . (1.1) 



For distances larger than r*, gravity is four-dimensional linearized Brans-Dicke, with pa- 
rameter uj = 0. Thus, a marked departure from conventional physics continues down to 
distances much smaller than ro, the distance at which the extra dimension is naively hidden. 
That departure is well-characterized and provides a possible signature for the existence of 
extra dimensions. 

There is, however, a catch. The cosmological solutions that drive interest in DGP gravity 
indicate that tq should be close to today's Hubble radius. Localized matter sources are 
embedded in a cosmological spacetime. Far enough away from a given source, the motion 
of observers and test particles is dominated by the cosmology, rather than the gravity of 
the matter source. However, if tq is today's Hubble radius, then the distance away from a 
source at which cosmology dominates the metric is also r^,. So, the DGP departure from 
Einstein gravity that was calculated in a static Minkowski background is only significant on 
length scales where the background cannot actually be approximated as Minkowski. The 
calculations needs to be refined if we are to have the correct new physics signature. 
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The subject of this paper is to look at corrections to Einstein gravity for the DGP 
braneworld theory in a cosmological background.0 We begin by laying out the necessary 
details of the model as well as the particular background of interest. We then solve for 
the metric of a spherically symmetric, static matter source in that background. We indeed 
find that the corrections to Einstein are sensitive to the background cosmology, even in the 
region inside where the cosmological flow is ostensibly irrelevant. We discuss possible 
tests for the detection of new physics at astronomical scales and suggest that it is possible 
for such tests of local gravity to reveal information about global features of the full five- 
dimensional cosmology, and ultimately shed light on the nature of dark energy and today's 
cosmic expansion. 



II. PRELIMINARIES 
A. The Model 



Consider a braneworld theory of gravity with an infinite-volume bulk and a metastable 
brane graviton We take a four- dimensional braneworld embedded in a five-dimensional 
Minkowski spacetime. The bulk is empty; all energy-momentum is isolated on the brane. 
The action is 

^(5) = - Y^^' / d'x^ R + J d^x^/^ Cm + Sgh • (2.1) 

The quantity M is the fundamental five-dimensional Planck scale. The first term in Eq. 
corresponds to the Einstein-Hilbert action in five dimensions for a five-dimensional metric 
qab (bulk metric) with Ricci scalar R. The term Sgh is the Gibbons-Hawking action. In 
addition, we consider an intrinsic curvature term which is generally induced by radiative 
corrections by the matter density on the brane [0]: 

-T^M| / rf^xy^ . (2.2) 

iOTT J 

Here, Mp is the observed four-dimensional Planck scale (see |1|,0J^ for details). Similarly, 
Eq. ( |2.2| ) is the Einstein-Hilbert action for the induced metric gj^J on the brane, R^^^ being 
its scalar curvature. The induced metric iM 



The work in this paper, as well as that in Refs. |jT2,13,rj], addresses issues similar to those studied 



in Refs. []lqjl7[| . Our results differ from those in the latter work because we are only concerned with 
solutions that asymptotically approach well-behaved, familiar solutions (e.g., Minkowski space or 
deSitter expansion) far away from the matter source, including far away from the source in the 
direction into the bulk. 

^Throughout this paper, we use A, B, . . . = {0, 1, 2, 3, 5} as bulk indices, i^, . . . = {0, 1, 2, 3} as 
brane spacetime indices, and i, j, . . . = {1, 2, 3} as brane spatial indices. 
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g^^ = d.X^d^X^gAB , (2.3) 
where X'^(x^) represents the coordinates of an event on the brane labeled by x^. The action 



given by Eqs. (|2.1|) and (|2.2|) leads to the following equations of motion 



1 ^ r/i \^('4) Svr 
—Gab + d{hia.ne)G\'B = —. 
2ro Mp 



Gab + 5(brane)G^]j = -—tTabI brane , (2.4) 



where G^b is the bulk Einstein tensor, G^^ is the intrinsic brane Einstein tensor, and 
^.Aslbrane IS the matter energy-momentum tensor on the brane, and we have defined a 
crossover scale 

^« - H ■ (2.5) 

This scale characterizes that distance over which metric fluctuations propagating on the 
brane dissipate into the bulk 



B. Field Equations 

We are interested in finding the metric for static, compact, spherical sources. We are 
interested in looking at these solutions in a cosmological background rather than a Minkowski 



background [|13|,|T4| to ascertain what affects cosmology might have on the observability of 
corrections to Einstein gravity. We restrict ourselves to a background deSitter cosmology. 
Not only is such a background the simplest and most convenient, but observations suggest 
that we are currently undergoing deSitter-like cosmic acceleration. If the Hubble scale of 
such an acceleration is varying slowly, the results obtained here would apply. They will 
also shed light on the technical VDVZ problem and how one recovers Einstein gravity in a 
deSitter background. 

Under this circumstance of a static spherical source in a deSitter background, one can 
choose a coordinate system in which the cosmological metric is static (i.e., has a timelike 
Killing vector) while still respecting the spherical symmetry of the matter source. Let the 
line element be 

ds^ = N\r, z)de - A\r, z)dr^ - B\r, z)[de^ + sin^ ed<P^] - dz^ . (2.6) 

This is the most general static metric with spherical symmetry on the brane. The bulk 
Einstein tensor for this metric is: 
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The prime denotes partial differentiation with respect to r, whereas the dot represents partial 
differentiation with respect to z. 

We wish to solve the five-dimensional field equations, Eq. ( p.4|) . This implies that all 
components of the Einstein tensor, Eqs. ( p.7|) , vanish in the bulk but satisfy the following 
modified boundary relationships on the brane. Fixing the residual gauge B\z=o = r, when 
z = and imposing 22^symmetry across the brane 
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(2.8) 



from Gtt, Grr, and Gqq, respectively. 
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C. Background Cosmology 



The deSitter solution with Hubble scale, H, has the following metric components:! 



N{r,z) = {lTH\z\){l-H^Py^^ (2.9) 
A{r,z) = {1 T H\z\) (l - H\^Y'^^ (2.10) 



B{r,z) = {It H\z\) r . (2.11) 
The brane energy-momentum tensor required for a given Hubble parameter is 

^ilbrane = S (z) diag (pH , -pH , -PH , -PH , 0) , (2.12) 

with 

PH = -PH = ^^iroH±l) . (2.13) 

For physical reasons, we restrict ourselves to both non-negative pn as well as non-negative 
H. Under these circumstances, one can see that H > r^^ in the self-accelerating phase. 
One can show that the solution Eqs. (|2.9| - |2.13D corresponds to a coordinate transformation 



of the deSitter solution in homogeneous cosmological coordinates found in Ref. 0]. 

III. SPACETIME GEOMETRY 

We have chosen a coordinate system, Eq. ( p. 61 ), in which a compact spherical matter 
source may have a static metric, yet still exist within a background cosmology that is non- 
trivial (i.e., deSitter expansion). Let us treat the matter distribution to be that required for 
the background cosmology, Eq. (|2.12| - ^rT3D , and add to that a compact spherically symmetric 



matter source, located on the brane around the origin {r = 0, z = 0) 

Ibrane = S{z) diag {pg{r) + pff, -Pg{r) + Ph, -Pgir) + Ph, -Pgir) + Ph, 0) , (3.1) 

where pg{r) is some given function of interest and Pg{r) is chosen to ensure the matter 
distribution and metric are static. 



The upper and lower signs in Eqs. (2.9-2.11) correspond to the two distinct cosmological phases 
that may exist for this theory Q . The upper sign corresponds to the choice of Friedmann-Lemaitre— 
Robertson-Walker (FLRW) cosmological phase, where a bulk observer views the braneworld as 
a relativistically expanding bubble from the interior; the lower sign corresponds to the self- 
accelerating cosmological phase, where the bulk observer views the braneworld as a relativistically 
expanding bubble from the exterior I^ IC]. These phases are vastly different geometrically and may 
have drastically different cosmological histories. We use this sign convention throughout the paper. 
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We are interested only in weak matter sources, Pg{r). Moreover, we are most interested 
in those parts of spacetime where deviations of the metric from Minkowski are small. Then, 
it is convenient to define the functions {n(r, z), a{r, z), b{r, z)} such that 



N{r, z) = 1 + n{r, z) 
y4(r, z) = 1 + a{r, z) 
B{r, z) = r [1 + b{r, z)] . 



(3.2) 
(3.3) 
(3.4) 



We may rewrite Eqs. ( p.7|) and (|2^ ) using the functions {n(r, z), a{r, z), b{r, z)}. 

In order to determine the metric on the brane, we will implement the approximation 



n 



z=0 



(3.5) 



even in the presence of a compact matter source. Equation ( |3.5| ) presumes that the contri- 
bution to n\z=o from the matter source is negligible compared to other terms in the brane 
boundary conditions, Eqs. (|2.8| ). With this one specification, a complete set of equations, 
represented by the brane boundary conditions Eqs. ( p.8| ) and Gl = 0, exists on the brane so 
that the metric functions may be solved on that surface without reference to the bulk. We 
justify this assumption by examining the full bulk metric in detail in the Appendix. 
The brane boundary conditions Eqs. (|2.8|) now take the form 



-(a + 26) = ro 



-2b = ro 



-{d + b) = ro 



2a 

~~2 



2a' 

2n' 2a 

n' 
r 



n" + 



+ ^i?;(r) + 3i/(ro//±l) 
//(3roi/ ± 2) 
+ Hi^'ir^H ± 2) , 



(3.6) 



where we have defined an effective Schwarzschild radius Rg{r) for the matter source inside 



a given distance, r. 



Rn 



dr r'^Pg{r) 



(3.7) 



and where we have neglected second-order contributions (including those from the pressure 
necessary to keep the source static). Covariant conservation of the source on the brane 
allows one to ascertain the source pressure, p(r), given the source density p(r): 



Pa 



-n pg 



(3.8) 



One can show that the Gz^^component of the bulk is identically zero on the brane when 
covariant conservation of the matter source and the brane boundary conditions Eqs. 
are satisfied. This observation is a restatement of one of the Bianchi identities. 
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From the brane boundary conditions Eqs. (|3.6|) , one may eliminate a and 6, and arrive 
at an equation which may be integrated imediately with respect to the r-coordinate. The 
result yields 



r n + ra H — r H 
2 



■g ' 



(3.9) 



where we have imposed a boundary condition requiring that the metric not be singular at 
the origin. Applying Eqs. ( |3.9| ) and (|3^ ) to the G^^-component of the bulk metric, one can 
again arrive at an equation that is dependent only on variables on the brane itself which is 
integrable with respect to the r-coordinate. Defining the quantity 



/(r) = rn' — 



{2roH ± 1) 



the integral of the equations 



2r 2ro 
on the brane yields 



(3.10) 



4r^ 
r 



— Rn 

2 ^ 



3^ 
2^0 







(3.11) 



where we have applied a second spatial boundary condition requiring that the metric not 
have spurious cusps at the origin. Equation ( |3.11| ) has two solutions for /(r). We select the 
solution that matches onto the proper background deSitter behavior for large-r. Then, 



fir) 



-(3 ± 2roH) + 3(1 ± 2roH) 



\ 



1 + 



9r3 (1 ± 2roHy 



(3.12) 



Substituting this back into Eq. (|3.10| ) and Eq. (|3^) , we may articulate expressions for n{r) 
and a(r). 



rn 



2r 



[1 + 5{r)] - H^r^ 



a = ^ [1 - 5(r)] + \h\^ 

where we have defined the quantity 6{r) such that 

3^3 



(3.13) 
(3.14) 



Sir) 



[I ± 2roH) 



\\ 9r3 (l±2roif)2 



(3.15) 



These expressions are valid on the brane when r <C r^, H~^. In both expressions, the first 
term represent the usual Schwarzschild contribution with a correction governed by 5(r) 
resulting from brane dynamics, whereas the second term represents the leading cosmological 
contribution. Let us try to understand the character of the corrections. 
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IV. GRAVITATIONAL REGIMES 



There are important asymptotic hmits of physical relevance for the metric on the brane 
Eqs. ( p.l3| - pTT^ ). We have two system parameters, the crossover scale, tq, and the cosmolog- 



ical horizon radius, H~^, and we wish to understand how spacetime geometry differs when 
each is much larger than the other, as well as when they are the same order of magnitude. 
When roH <C 1, we recover familiar results for a compact source in a Minkowski back- 



ground [13,n|. The deSitter background becomes irrelevant at scales r <^ tq. A character- 
istic distance {rgRgY^^ demarks an Einstein phase close to the source from a Brans-Dicke 
{uj = 0) phase farther away. 



A. Einstein Phase 



Cosmological effects become important to the metric Eqs. ( p.l3| - pTT^ ) when roH ~ 1 or 



tqH ^ 1. Typically, we are concerned with the details of the metric when the influence of 
a given compact matter source dominates the local geometry. The competition between the 
leading Schwarzschild term, ~ Rg/r, versus the leading cosmological contribution, ~ H^r^, 
implies that when 

/R \^/^ 

r « (^) , (4.1) 

the local source dominates the metric over the contributions from the cosmological flow. In 
this region, Eqs. ( p.l3| - |3rT^ ) reduce to 




n = -^±Ji (4.2) 

(4.3) 

Notice that, indeed, there is no explicit dependence on the parameter governing cosmological 
expansion, H. However, the sign of the correction to the Schwarzschild solution is depen- 
dent on the global properties of the cosmological phase. The FLRW phase has a potential 
steeper than Einstein gravity, whereas the self-accelerating phase has a potential shallower 
than Einstein gravity. Thus, we may ascertain information about bulk, five- dimensional 
cosmological behavior from testing details of the metric where naively one would not expect 
cosmological contributions to be important. Indeed, for example, one can use local gravity 
tests to distinguish whether some local brane vacuum energy or self-acceleration is the cause 
of today's cosmic accelerated expansion. 
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B. Weak Brane Phase 



Even when the cosmological flow dominates the metric, one can still examine the per- 
turbed effect a matter source has in this region. When, r > {RgH-'^f/^, while still m a 
region well within the cosmological horizon (r ^ H~^), 



n 



Rg_ 

2r 



1 + 



1 



Rg^ 

2r 



3(1 ± 2roH) 
1 

3{l±2roH) 



2 



(4.4) 
(4.5) 



This is the direct analog of the weak-brane phase one finds for compact sources in Minkowski 
space. The residual effect of the matter source is a linearized scalar-tensor gravity with 
Brans-Dicke parameter 



(4.6) 



Notice that as TqH —>■ oo, we recover the Einstein solution, corroborating results found for 
linearized cosmological perturbations 



mi 



Moreover, note that in the self-accelerating 
cosmological phase, the scalar component couples repulsively, though recall from Eq. (|2.13|) 
that H > Tq ^ in this phase. 



C. The Picture 



One can consolidate these results and show from Eqs. ( p.lSHOSi ), that there exists a 
scale, 



'l±2r^Hf 



1/3 



(4.7) 



inside of which the metric is dominated by Einstein but has corrections which depend on the 
global cosmological phase, i.e., Eqs. (|4.2|-|i7^). Outside this radius (but at distances much 
smaller than both the crossover scale, tq, and the cosmological horizon, H~^) the metric is 
weak-brane and resembles a scalar-tensor gravity in the background of a deSitter expansion, 
i.e., Eqs. (OHJ). 

The picture one develops is that the metric is scalar-tensor, but with a radially dependent 
coupling of the scalar to matter. As one gets closer to a source (as gravity gets stronger), 
the scalar coupling gets suppressed. There are two basic tests that we can pose. First, small 
corrections to the Newtonian potential using Eq. ([4.2|) and second, comparison of the ratio 
of the gravitomagnetic to gravitoelectric (Newtonian) forces. This latter may be encoded 
by the effective Brans-Dicke parameter, Eq. 
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V. PHYSICAL CONSIDERATIONS 



A. Constraints on tq 

It is useful to review the constraints on the crossover scale, tq. The constraint may be 
written as a constraint on the fundamental Planck scale: 

10"=^ eV < M < 1 GeV . (5.1) 

If M is too small, quantum gravity effects become important at short distances. The lower 
bound is provided by constraints from millimeter tests of Newtonian gravity and other 
constraints coming from loss of energy into the extra dimension . Cosmology provides the 
upper bound for M. If M is too large, the crossover scale becomes too small to account for 
the relationship between the observed Hubble scale and the independently measured matter 
density These bounds on the fundamental Planck scale correspond to 

1 Gpc < ro < 10^" Gpc . (5.2) 

An even more stringent case applies to the self-accelerating cosmological phase. In this 
phase, the Hubble scale is bounded from below H{t) > r^^, where at late times, the uni- 
verse is deSitter and H r^^. If one goes further and postulates that the current cosmic 
acceleration is caused entirely by this late-time self-acceleration, then, using constraints from 
Type lA supernovae 0, the best fit for tq is 

ro = 1.2ltH^^o-' , (5.3) 
where Hq is today's Hubble scale. Taking Hq ^ 70 km s~^Mpc~^, 

ro ~ 5 Gpc . (5.4) 
We are particularly interested this possibility and, in this Section, take ro to have this value. 

B. Gravitational Lensing 

The lensing of light by a compact matter source with metric Eq. ( ^.13| - P?T^ ) may be 
computed in the usual way. The angle of deflection of a massless test particle is given by 

= [dr 4 ; ^ , (5.5) 

where E = N'^dt/dX and J = r'^d(f)/d\ are constants of motion resulting from the isometrics, 
and d\ is the differential affine parameter. Then for any metric respecting the condition 
Eq. ( |3.9|) , the angle of deflection is 
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FIG. 1. Mass discrepancy, 5M, for a static point source whose Schwarzschild radius is r„. The 



sohd curve is for a self-accelerating background with H 
background with H = r^"^. 



The dashed curve is for a FLRW 



TT 
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(5.6) 



where 6 is the impact parameter. This result is equivalent to Einstein, so we see that 
light deflection is unaffected by DGP corrections. This is consistent with the idea the 
DGP corrections correspond to an anomalous spatially-dependent scalar coupling. Since 
scalars do not couple to light, the trajectory of light in a gravitational field should remain 
unaffected. Then, lensing measurements probe the true mass of a given matter distribution. 
One can then compare that mass to the mass taken from assuming the gravitational force 
is Newtonian. 

The mass discrepancy between the lensing mass (the actual mass) and that determined 
from the Newtonian force may be read directly from Eqs. (p.l3| ) and ( p.l5| ). 

Miens ^ 1 



5M 



1 



1 . 



(5.7) 



MNewt ' 1 + 5(r) 

This ratio is depicted in Fig. for both cosmological phases. When the mass is measured 
deep within the Einstein regime, the mass discrepancy simplifies to 

/ .3 \ 1/3 



5M 



(5.8) 



Solar system measurements are too coarse to be able to resolve the DGP discrepancy between 
lensing mass of the sun and its Newtonian mass. The discrepancy 6M for the sun at 0{AU) 
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scale distances is approximately 10~^^. Limits on this discrepancy for the solar system as 
characterized by the post-Newtonian parameter, 7 — 1, are only constrained to be < 3 x 10"'^. 

A possibly more promising regime may be found in galaxy clusters. For 10^^ — > 10^^ Mq 
clusters, the scale (rQRgY^^ has the range 6 — 14 Mpc. For masses measured at the cluster 
virial radii of roughly 1^3 Mpc, this implies mass discrepancies of and taking their virial 
radii to be roughly 5 — > 8%. X-ray or Sunyaev-Zeldovich (SZ) measurements are poised to 
map the Newtonian potential of the galaxy clusters, whereas weak lensing measurements 
can directly measure the cluster mass profile. Unfortunately, these measurements are far 
from achieving the desired precisions. If one can extend mass measurements to distances on 
the order of tq, Fig. |^ suggests discrepancies can be as large as —10% for the FLRW phase 
or even 50% for the self-accelerating phase. 



C. Orbit Precession 



Imagine a body orbiting a mass source where Rg{r) = Tg = constant. The perihelion 
precession per orbit may be determined in the usual way 

f , J AN , , 

M= jdr - ^ ^ ^ , (5.9) 



where E = N'^dt/ds and J = r'^dip/ds are again constants of motion resulting from the 
isometrics, and now ds is the differential proper time of the orbiting body. Assuming a 
nearly circular orbit and that we are deep within the Einstein regime (so that we may use 
Eqs. (13-13)), then 

A^ = 2, + ^T^f^)'". (5.10) 

The second term is the usual Einstein precession. The last term is the new anomalous 
precession due to DGP brane effects. Note that this correction is the same as one would get 
if one assumed a purely Newtonian potential Eq. ( ^.2| ) without spatial metric effects. The 
correction to the precession rate one expects from DGP gravity is 



— = T — 
dt 8ro 



DGP - =F— = =F5 yuas/year . (5.11) 



Note that this result is independent of the source mass, implying that this precession rate 
is a universal quantity dependent only on the graviton's effective linewidth (vq^) and the 
overall cosmological phase. Compare Eq. ( |5.11D to the classic Einstein precession correction 
for nearly circular orbits: 
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For increasing r, the distance from the sun at which the DGP correction begins to overtake 
the first Einstein correction is 37 AU. 

Nordtvedt |^ quotes precision for perihehon precession at 430 /las/year for Mercury 
and 10 yuas/year for Mars. Improvements in lunar ranging measurements suggest 
that the Moon will be sensitive to the DGP correction Eq. ( ^.11 ) in the near future. Also, 



BepiColombo, an ESA satellite being sent to Mercury at the end of the decade, will also 
be sensitive to this correction ||2^. Incidentally, it is interesting to contrast these numbers 
with a precision of 4 x 10^ /ias/year for the rate of periastron advance in Binary Pulsar 
PSR 1913+16 p3. The solar system seems to provide the most promising means to constrain 



this anomalous precession from DGP gravity. 

VI. CONCLUDING REMARKS 

The braneworld theory of Dvali-Gabadadze-Porrati (DGP) is an intriguing extension 
of Einstein gravity that exploits the possible existence of infinite-volume, extra dimensions. 
It is a theory where the four-dimensional graviton is effectively metastable, and provides 
a novel alternative to conventional explanations of the dark energy that is responsible for 
today's cosmic acceleration. 

In this paper, we detailed the solution of static, spherical matter sources in the back- 
ground of deSitter cosmology for DGP gravity. The gravitational field of a matter source 
exhibits important dependences on cosmology. Residual dependences on the full five- 
dimensional cosmological phase also exist in the regime deep in the gravity well of the 
matter source where the effects of cosmology are ostensibly irrelevant. These residual de- 
pendences allow one to use local (e.g., solar system) measurements of the gravitational field 
to ascertain details of the global cosmology. 

In DGP gravity, we find that massless test particles can probe the true mass of a matter 
source, whereas tests of the source's Newtonian force leads to discrepancies with general 
relativity. These discrepancies translate into a universal anomalous precession, as large 
as ±5 /ias/year, suffered by all orbiting bodies. The numerical value of this anomalous 
precession is dependent only on the graviton's effective linewidth and the global geometry 
of the five-dimensional cosmology. Current constraints on Mars' orbit are on the threshold 
of being sensitive to this anomaly . Future improvements in lunar ranging |^ as well 



as data from satellite missions at the end of the decade p3[ should be sensitive to possible 



corrections due to DGP and gravitational leakage into extra dimensions. 
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APPENDIX 

In order to see why Eq. ( p.5| ) is a reasonable approximation, we need to explore the full 
solution to the bulk Einstein equations, 

GAB{r,z) = Q, (Al) 



satisfying the brane boundary conditions, Eqs. ( p.8| ), as well as specifying that the metric 
approach the deSitter solution Eqs. ( |2.9| - |2.1I| ) for large values of r and z, i.e., far away from 
the compact matter source. 

First, it is convenient to consider not only the components of the Einstein tensor 
Eqs. ( [2. 71 ), but also the following components of the bulk Ricci tensor (which also vanishes 
in the bulk): 



'N" N' A' ^N' B' 



+ 



N NA N B 

— + + 2 

N N A N B 



(A2) 



N A 2B , 

«' = iv + i + ir- 

We wish to take G^r = 0, = 0) ^"^^ — ^ derive expressions for A{r, z) and -B(r, z) 
in terms of N{r,z). Only two of these three equations are independent, but it is useful to 
use all three to ascertain the desired expressions. 

Since we are only interested in metric when r, z <^ tq, H"^ for a weak matter source, 
we may rewrite the necessary field equations using the expressions Eqs. ( p.3|) . Since the 
functions, {n(r, z), a{r, z), b{r, z)} are small, we need only keep nonlinear terms that include 
z-derivatives. The brane boundary conditions, Eqs. ( p.8| ), suggest that a and b terms may 
be sufficiently large to warrant inclusion of their subleading contributions. It is these z- 
derivative nonlinear terms that are crucial to the recover of Einstein gravity near the matter 
source. If one neglected these bilinear terms as well, one would revert to the linearized, 
weak-brane solution (cf. Ref. [0). 

Integrating Eq. ([A3|) twice with respect to the 2;-coordinate, we get 
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n + a + 2b = zgi{r) + g2{r) , (A4) 



where gi{r) and g2{r) are to be specified by the brane boundary conditions, Eqs. ( p.8|) , and 
the residual gauge freedom 6b{r)\z=o = 0, respectively. Integrating the G^^-component of 
the bulk Einstein tensor Eqs. (|2.7|) with respect to the z-coordinate yields 



r{n + 26)' - 2 (a - 6) = g-iir) . (A5) 

The functions gi{r), g2{i^), and g3{r) are not all independent, and one can ascertain their 
relationship with one another by substituting Eqs. ( [A4D and ( [A5D into the bulk equation. 
If one can approximate n = for all z, then one can see that Gzr = 0, = 0, and = 
are all consistently satisfied by Eqs. ( |A^ ) and (|A5|) , where the functions gi{r), g2{r), and 
(73 (r) are determined at the brane using Eqs. ( |3.13|) and ( |3.14] ) and the residual gauge freedom 
%)U=o = 0: 



g,{r) = T^H - {R,6y (A6) 
92{r) 



^'n-6)+ r dr ^{1 + 5) (A7) 
JO 

g,{r) = ^{l - ?,6) - 2Hh\ (A8) 
2r 



2r 



where we have used the function 5(r), defined in Eq. ( p.l5| ). Using Eqs. (|A4|-|A8|) , we now 



have expressions for a(r, z) and 6(r, z) completely in terms of n(r, z) for all (r, z). 

Now we wish to find n{r, z) and to confirm that h = is a good approximation 
everywhere of interest. Equation ( [A^ ) becomes 



2??' 

n" + — + fi = ±H[gi (r) ± if] , (A9) 

where again we have neglected contributions if we are only concerned with r, z <^ tq, H~^. 
Using the expression Eq. ( |A6|) , we find 

n" + ^ + n = -3H' t ^ [i?,<5(r)]' . (AlO) 

Then, if we let 

1 r I 

n = lTHz ifV^ =F roH / dr —RJr)S(r) + 5n(r, z) , (All) 

2 JO r"' 

where 5n(r, z) satisfies the equation 

Sn" + ^ + 5n = , (A12) 
r 

we can solve Eq. (|A12|) by requiring that 5n vanish as r, z ^ cxd and applying the condition 
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r 5n'l=o = ^ [1 + (1 ± 2r,H)5{r)] , (A13) 

on the brane as an alternative to the appropriate brane boundary condition for 6n{r, z) 
coming from a hnear combination of Eqs. ( p.8|) . We can write the solution explicitly: 

/■oo 

5n{r,z)= f dk c{k)e~^'' sinkr , (A14) 
Jo 



where 



2 

c{k) = — I dr rsinkr 5n|^^Q (r) . (A15) 



vr Jo 

We can then compute 5n\z=o, arriving at the bound 



r Jo 

for all r -C ro, H~^. Then, 



SnU=o<z [dr^, (A16) 



h\,=o = TH + 5n\,=o ■ (A17) 



When the first term in Eq. ( [A 17]) is much larger than the second, Eq. (|3.5|) is a good 



approximation. When the two terms in Eq. ([A17| ) are comparable or when the second term 



is much larger than the first, neither term is important in the determination of Eqs. (|3.13|) 



and ( |3.14| ). Thus, Eq. ( p.5|) is still a safe approximation. 

One can confirm that all the components of the five- dimensional Einstein tensor, 
Eqs. ( p.7|) , vanish in the bulk using field variables satisfying the relationships Eqs. ( [A4D , 



and ( [All| ). The field variables a{r,z) and b{r,z) both have terms that grow with 



z, stemming from the presence of the matter source. However, one can see that with the 
following redefinition of coordinates: 

R = r- zro^ (A18) 

Z = z+rdr^, (A19) 
JO 

that to leading order as 2; H^^, the desired Z-dependence is recovered for a{R, Z) and 
b{R, Z) (i.e., ^HZ), and the Newtonian potential takes the form 

n{R, Z) = tHZ - ^H^R^ + ■■■ . (A20) 

Thus, we recover the desired asymptotic form for the metric of a static, compact matter 
source in the background of a deSitter expansion. 
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